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UJ ' We determine the full mass and q^ dependence of the heavy quark vacuum polarization function 

^■^ • n(g^) and its contribution to the total e^e~ cross section at 0{a^) and 0{al) in perturbative 

CN| ' QCD. We use known results for the expansions of n(q^) at high energies, in the threshold region 

and around q^ = 0, conformal mapping and the Pade approximation method. From our results for 

n(g^) we determine numerically at 0{al) the previously unknown non-logarithmic contributions in 

f~ i . the high-energy expansion at order (m^/g^)' for i — 0,1 and the coefficients in the expansion around 

Mh' g^ = at order g^" with n > 2. We also determine at 0{a^) the previously unknown 0{v") constant 

C») ; term in the expansion of Tl{q^) in the threshold region, where v is the quark velocity. Our method 

Q^ ' allows for a quantitative estimate of uncertainties and can be systematically improved once more 

i-i^ . information in the three kinematic regions becomes available by future multi-loop computations. 

For the contributions to the total e+e" cross section at 0{a'l) we confirm results obtained earlier 

^S) ■ by Chetyrkin, Kiihn and Steinhauser. 
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I. INTRODUCTION 



The vacuum polarization function n(g^) defined by the correlator of two electromagnetic currents j^{x) 
-0(a;)7'V(a;), 



(.9^1''?^ - Qfj-Qi^] 



n(q2) = -^Jdxe^'^-{0\TJ^{x)JAO)\0) , (1) 



where q^ is the four- momentum of the quark pair produced or annihilated by j^, represents an important quantity 
for theoretical studies as well as for many practical phenomenological applications. Relevant applications for the case 
of massive quarks include predictions of the hadronic cross section R ~ Im[n], or sum rules for the determination 
of the heavy quark masses |l|, y]. These sum rules are based on moments of the cross section for heavy quark pair 
production 

f°° ds 

Mn^ ^^Ris), (2) 

which in fixed-order perturbation theory are related to the expansion coefficients of Tl{q^) around q^ = 0, 

_j oo 

n('Z^-0,m^) = — — ^^M„g^", (3) 

where Qq is the heavy quark electric charge. In general the knowledge of the full dependence of the vacuum polarization 
function Tl{q^) on q^ and the quark mass m is desirable to avoid having to rely on approximations that are only valid 
in certain kinematic regimes. 

At 0{as) the full mass and q^ dependence of the vacuum polarization function is known from analytic computations 
carried out in Rcf. Isl. At 0{a^) analogous analytic results exist for the contributions that originate from inserting 
the massive [0, 0] and massless pi fermion loops into ©(a,,) onc-gluon exchange diagrams. For the other 0{a^) 
contributions results for the expansions of 11(9^) in the high-energy limit, \q'^\ — > oo, the nonrclativistic threshold 
regime, q^ w 477i^, and in the Euclidean region around q^ = were used to reconstruct an accurate approximation 0, ^. 
The method is based on the definition of subtraction functions which account for all logarithmic terms that arise for 
the expansions in the high-energy limit and in the threshold region. Using a conformal transformation to a new 
variable to the full q^ and mass dependence in the complex plane of the remaining contributions can be mapped 
into the unit circle rendering those contributions to an analytic function in the variable uj. The latter can then be 
successfully approximated by Fade approximants using the remaining expansion coefficients that are not related to 
logarithmic terms. With a large number of expansion coefficients for the three kinematic limits the full mass and g^ 
dependence of the 0{al) vacuum polarization function can be determined with small numerical uncertainties. 

For the 0{al) vacuum polarization there is also no fully analytic result available in the literature. A numerical 
study of the full 0{aln'i) double fcrmionic contributions to the vacuum polarization function can be found in Rcf. [Bl. 
In the high-energy expansion its contributions to the total cross section up to order (?7i^/g^)^ are known [Q. A 
comprehensive review of these results can be found in Ref. |ll[. Moreover, in the threshold region, where an expansion 
in the small quark velocity v can be carried out, the 0{al) contributions to the total cross section at order 1/w' for 
i = 2, 1,0 are available from a factorization theorem for the heavy quark-antiquark pair production cross section in 
nonrclativistic QCD (NRQCD) at next-to- next-to-leading order (NNLO) [|2|, |3|. More recently also the moments 
Ml llJ, |l^ and M2 Q at 0(0;^) have become available using elaborate high-power computer algebra tools. 

In this work wc use the presently available information on the C(a^) corrections to the vacuum polarization function 
n(g^) in the high-energy limit, in the threshold region and the small q^ domain to reconstruct the full q^ and mass 
dependence of the vacuum polarization function at 0{al). The method we use is similar to the approach of Rcfs. [^ ||] 
employed previously for the C(ag) corrections of the vacuum polarization function (see also Rcfs. [|l^ 18, |l^), but 



also accommodates a few notable differences which are motivated by the fact that less information is known on the 
vacuum polarization function at 0{al). While in Refs. [^ || information from the high-energy expansion up to order 
{w? /q^) and from the threshold expansion up to next-to-leading order (NLO) was incorporated for the reconstruction, 
we account for the expansions up to order {m? /q^Y ^^ high energies and up to NNLO in the threshold region. While 
in Ref. H, p| the full set of terms in the high-energy expansion of 11 was used for the construction, in this work we 
only rely on the terms that carry an absorptive part and that contribute to the cross section above threshold. We 
show that our method allows to determine previously unknown non-logarithmic terms of the vacuum polarization 
at 0{a'l) in the high-energy expansion with very small errors. Moreover, while in Ref. [M 0] and later in Ref. pH] 



the coefBcients of the small-g^ expansion were included up to order g^** and q^^, respectively, we only rely on the 
presently available 0{al) coefficients up to order q^. Our method allows to determine the expansion coefficients at 
order g^" with n > 3. The results allow to compute the corresponding moments M„ in the fixed-order expansion at 
0{a^). For phenomenologically relevant values of n the error in the M„ due to the uncertainties in these coefficients 
is an order of magnitude smaller than the remaining scale- uncertainties of the M„ at 0{al) . We demonstrate the 
reliability of the results by using the same approach for determining the corresponding coefficients for the vacuum 
polarization function at 0(0^.) where their values are well known analytically from the computations of Feynman 
diagrams. Another noteworthy difference of our approach to Refs. H, M is that we implement a continuous set of 
subtraction functions to have a more reliable estimation of the uncertainty inherent to the method. Our approach 
can systematically incorporate new information from the expansions in the three kincmatical regions, once it becomes 
available. 

One important application of the vacuum polarization function at 0{a^) obtained in this work is an analysis of 
low-n moments of the e'^e~ —> cc cross section to determine the MS charm quark mass rUc and to investigate the 
uncertainty in rric that arises from the difference of using fixed-order and contour-improved perturbation theory. 
For using contour-improved perturbation theory, which involves integrations of n(g^) in the complex g^-plane, it is 
essential to have the full mass and q^ dependence of the vacuum polarization function. Such an analysis was carried 
out in Ref. ||2l| at 0{a^). Determinations of the charm quark mass rfic using the vacuum polarization function at 
0{ag) in the fixed-order expansion alone were carried out recently in Refs. [ |l5| , |22| . In this paper we discuss in 
detail the reconstruction of the full q^ and mass dependence of the 0{al) and 0{al) corrections to the vacuum 
polarization function as outlined above. The thorough analysis of uncertainties in the charm and bottom quark MS 
masses obtained from low-n moments of the e+e~ cross section will be given in a subsequent publication. 

The program of this paper is as follows: In Sec. O we set up our notation and in Sec. Ill wc present the basic 



features of our method for reconstructing the vacuum polarization function. In Sec. IV we explain details about 
how logarithmic contributions in the expansions in the threshold region and for high energies are incorporated and 
in Sec. M we present how the remaining non-logarithmic terms are treated. Some of the solutions we obtain have 



unphysical properties. Criteria that allow to identify and discard such solutions are discussed in Sec. VI. Numerical 
analyses for the 0{al) and 0{al) contributions of the vacuum polarization function are given in Sees. VII and VIII. 
Our conclusions are given in Sec. px| . 

II. NOTATION 

The relation between the normalized e~^e~ cross section R and and the vacuum polarization function 11 reads 

i?(g2) = 127rg2lmn(q2+i0,m2), (4) 

where Qq is the heavy quark electric charge. The perturbative fixed-order expansion of Il{q^,Tn^) has the form 



n(g2,m2) =n(")(g2,m2) + f££^Al^\ Il^^\q^,m^ 



"^('^ )^ ^(^)('Z^-^^^) + f^f^] U^^Hq^m^^') + ... , (5) 



with the color factor Cp — 4/3. We use the on-shell normalization of the vacuum polarization function, where 

n(0,TO2) ^0. (6) 

We exclude the so-called singlet contributions where the vacuum polarization function contains a three-gluon cut. 
Note that in this work we do not distinguish between the contributions in II'^' and II^'^^ proportional to the different 
SU(3) group theory color factors since there isn't any compelling technical reason that would make such a distinction 
mandatory. This approach neglects the existence of the multi-particle cuts from diagrams with the insertion of 
massive fermion loops. Their contribution is strongly phasc-space-supprcssed and can be safely ignored for the level 
of accuracy intended in this work. We emphasize, however, that our approach can be applied to the individual color 
contributions as well. 

For the reconstruction of the vacuum polarization function accomplished in this work we use exclusively the pole 
mass scheme, m = mpoic, since it allows for the most transparent treatment of the information from the quark pair 



production threshold. Moreover we use the choice ji ^ m = rripoie for the renormalization scale and generally suppress 
the /i-dependence of the functions II'-'-'. To simplify the presentation we frequently use the variable 

(7) 



Am? 



For the strong coupling we use Uf ~ ni-\-\ active running flavors, where quarks that are heavier than those produced 
by the current j^ are integrated out and where all n^ light flavors are treated as massless. 

The analytic expression for the vacuum polarization functions at 0{as) [^ is an important ingredient of our analysis. 
The corresponding contributions using the notation of Eq. (^ have the form 



n(") = 


3 

167r2 


20 4 

y + 3] 


n(i) - 


3 

167r2 


■5 13 
6"^ 6z 



4(l-z)(l + 2z) 
3z 

(l-z)(3 + 2z) 



G{z) 



G{z) 



(l-z)(l-16z)^2 



6z 



G\z) 



(l + 2z) 
6z 



l + 2z(l-z) 



dz 



I{z) 



(8) 



where 



I{z) = 6 C3 + 4Li3(-u) + 2Li3(M) 
2 u Inu 



2Li2(-u) + Li2(u) 



lnM-2 



2 ln(l 



ln(l - u) 



In^' 



G{z) 



1 



with u 






(9) 



An important application is the determination of the moments Mn in the fixed-order expansion. Here the pole 
mass scheme is strongly disfavored since it contains an ©(Aqcd) renormalon ambiguity that leads to a quite bad 
perturbative expansion of the moments. For small values of n this problem can be avoided conveniently by using 
the MS mass scheme. For the complications that arise for large values of n see e.g. Refs. |^, |2J]. In this paper we 
use the MS running mass m with n^ = n^ + 1 running flavors for discussions of the moments A/„. Using a common 
renormalization scale /i for the mass and the strong coupling, the fixed-order perturbative expansion of the moments 
Mn can be written in the form [l„i^ = \n{m'^ (jj.) / jj^)] 



M„ 



Q: 



9 



4 (4TO2(/i))" 



■^(0) + ^ ((.aO) + ^(11)^^^) ^ (^^y (^(.0) ^ ^(21)^^^ ^ c!,22,Q 



sif^) 



f/=.(30) , /=.(31), , /=,(32),2 ,/j{33)j3 \ 



(10) 



adopting the notation of Refs. |j, ||. 



III. THE METHOD 



The expansions of Ii{z) in the threshold region z ~ 1 and the high-energy limit |z| — > cx) involve powers of log(l — z) 
and log(— 4z), respectively. Above production threshold for z > 1 these logarithmic terms contribute to the absorptive 
parts in 11 that constitute the cross section according to Eq. (^. On the other hand, the expansion around z = 0, 
which is located in the Euclidean region, leads to fully analytic terms and admits a usual Taylor expansion. We want 
to reconstruct the full q"^ and mass dependence of H^^^ by building functions that incorporate all known properties 
of n'^^^ in the threshold regime, the high-energy limit and the region around z = 0. We carry out the same program 
also for n^^^ using only coefficients in the expansions that are analogous to the available information for 11^'^^ From 
the reconstructed n^'^^ we can determine previously unknown non-logarithmic coefficients in the high-energy and the 
nonrelativistic expansions as well as the 0{a^) corrections of the moments M„ for n > 3. Using the reconstructed 
H?' function we can test the reliability of these determinations and find that these coefficients and moments can be 
determined remarkably well. 

Following the approach of Ref. M, we split Il?'^^{z) into two parts. 



n(2,3)( 



n(2.3) 



(^) 



t(2,3)^ 



n^r^z) 



(11) 



where IIj^' ' [z) are designed such that they contain the logarithmic terms in the expansions around z = 1 and for 
\z\ — > oo. They can be conveniently constructed from the functions II'^-'^) and G{z) given in Eqs. (H) and (O) since 
the latter readily incorporate analytic structures that allow to incorporate the appropriate threshold and high-energy 
behavior into nj^' \z). Once IIj ' (z) has been specified, the remaining task is to construct a Pade approximant for 

Ilrog that allows to incorporate the remaining non- logarithmic constraints in the regions z ~ 1, |z| -^ oo and z ~ 0. 
The general structure of a Pade approximant P„_m has the form 

En i 
r ^x) = '=0 °"^ fl9) 

1 + Ej=i bjx^ 

which means that there are n + m + 1 coefficients that need to be specified. Note that the coefficients Ui and bj are 

(2 3) • 

real numbers. Since Ilicg still has a physical cut for z > 1 along the positive real z axis, one cannot use the variable 
z to formulate the Pade approximant. A convenient variable to automatically account for this cut is to defined by (see 
e.g. Refs. 0, |l|) 



oj = , z = -77. (13) 

Here, the physical z-plane is mapped into the unit-circle of the complex w-plane, where approaching the physical cut 
from the upper (lower) complex z-half-plane corresponds to approaching the upper (lower) semi unit-circle in the 
complex tj-plane. The three points z = (0, 1, ±00) are conformally mapped onto uj = (0, 1, —1). Expressed in terms of 
the variable uj, Ilrcg can therefore be approximated by rational functions involving the Pade approximant P„_,„(ti>). 
All Pade approximants that turn out to have unphysical poles inside the unit circle have to be discarded. In practice 
some additional restrictive criteria have to be imposed to avoid an unphysical behavior of 11 and R due to poles in the 
Pade approximant outside the unit circle that are either close to the unit circle or have a large residue. We discuss 



these restrictions in Sec. VI 



(2^3)/ \ , ,, T3„jx n„^^^^A^n„+ f^v. tt(2,3) 



It goes without saying that the constructions of n[^^ (z) and the Pade approximant for Ilrog are not unique and 

that the resulting reconstructed fl'^'"^) functions have a dependence on choices made for their construction. The 

ambiguity in the procedure therefore needs to be quantified by accounting for variations in the construction. While 

in Ref. W p| variations coming from different choices for Pn,m were included for the error estimate, we include in 

(23) 
our work in addition continuous variations in the construction of IIj^' . We test the reliability of the method by 

determining properties of II'^' that are precisely known analytically, but that have not been incorporated for the 
construction of the approximation for 11^^^ 



IV. DESIGNING n,. 



To determine IIj^^ we need to account for the logarithmic terms that arise in n'^'^' in the threshold region 



z 



and in the high-energy limit |z| ^ cx). To facilitate the presentation it is convenient to write 

<'' (-) = ^it'' i^) + ^f i^) + ni^er^o' iz) , (14) 

where njj^j. ' and Il\^'i are designed to account for the logarithmic terms at threshold and at high energies, 

(2 3)' 

respectively, and Ilzcro incorporates subtractions that ensure a physical behavior at z = 0. 



Threshold Logarithms. We start by presenting the expansions of n(^^^'^^(z) and G{z) in the threshold hmit z ^ 1 
keeping terms up to NNLO in the expansion in \/l — z : 

n(i'(z) = - 0.1875 hr(l-z)- 0.314871 + 0.477465\/r^^ 

+ (0.354325 + 0.125 ln(l - z)) (1 - z) + O f (1 - zfl"^ 



1 700^7 

n(2) (z) = ' + (0.34375 - 0.0208333ni?) In2(l - z) + (0.0116822n<? + 1.64058) ln(l - z) + K^^'> 

\J\~ z 



0.721213 -0.0972614n£ + 3.05433 ln(l-z)\/T^^ + 0\{\-z) 



n(3)(z) 



2.63641 0.678207n(- - 27.2677 ,„ ,^, „ „ ,^, ,,logl-z 

— + -^ + (0.57419nf - 9.47414 '^ ) ' 

+ (-0.00231481n? + 0.0763889n<> - 0.630208) log3(l - z) 
+ (0.00194703nf + 0.0312341?7,^ + 1.3171) log2(l - z) 

+ (-0.06908487i| + 2.37068nf - 17.6668) log(l - z) + K^^^ + O ((1 - z)^^^) 



G{z) 



nVT 



2vr 



+ 0[{l-z)) 



(15) 



To avoid cluttering we show the various coefficients for Il^^''^'^'{z) in mmrerical form, but keep the number ni of light 
quark flavors as a variable. The expansions of II^^^ and G arc known from their exact expressions given in Eqs. (|8|) and 
(^ while the expansion for II'^^ can be derived from the results for R in the threshold region computed in Ref. P5[ . 
The expansion for II'^^ is obtained from the NNLO threshold cross section factorization formula for R within NRQCD 
first derived in Rcf. [|l^, Q (see also Ref. [^). The result was later confirmed by many other groups |g^. Note 
that within NRQCD it is the standard convention that only the rii light quark species contribute to the running of 
the strong coupling. Switching to nj = n£ + 1 running flavors affects the coefficient of the term ex ln(l — z) in H'-^K 



All other coefficients shown in Eq. (15) arc unaffected. We also note that the singlet contributions to the vacuum 
polarization function only affect the threshold expansion at N^LO in the expansion in yl — 2 and do not contribute 
at the order we consider here.^ The constant terms A'^^'"^^ that appea r in the nonre lativistic expansion of n^^'"^^ have 
not yet been computed from Fcynman diagrams. As wc show in Sees. VII and VIII they can be determined from the 
reconstructed vacuum polarization function based on the method described in Sec. III. 



ln™(l-z) 


[n(i)(^)]- 


{l-zy/^ln"'{l-z) 


[G{z)r[nW(z)r 


(l-z)"/2in-(l-z) 


{I - zr[G{zr[n^'\z)r 



TABLE I: First column: Logarithmic terms that arise in the expansion of II' ' (z) close to threshold where z ~ I. Second 
column: Corresponding functions used in the construction of H^j^^ (z)- 



(23) 
To construct II^j^j. we have to find appropriate functions that account for the different combinations of the loga- 
rithmic term ln(l — z) and powers of \/l — z that appear in Eqs. (p^. A convenient choice is given in Tab. |l and 
leads to 



t(2) 



^ItU^) - A^o 



(2) 



t(3) 



^IZi-) = K 



(3) 



1 + a^n^ z 



1 + ai^^ z 



n(i) (z)l + Af )n(i) (z) + 4^' (1 - z)G(z)n(i) (z) 



n(i)(z) 



A 



(3) 



1 + ai z 



n(i) (z) + Af i±^^^G(z)n(i) (z) + 4^^n(i) (z) , 



(16) 



^ Within NRQCD the dominant effect of the singlet contributions is associated to 4-quark operators with a Wilson coefficient that 
incorporates the hard effects of the 3-gluon annihilation. This operator leads to a momentum space potential oc a^/m'^. 



The coefficients A^ ' ' can be unambiguously determined from the expressions shown in Eqs. (p^). Obviously the 
choices in Tab. | are not unique. To have some quantitative way to account for this source of uncertainty we have 
multiplied the functions related to the highest power of ln(l — z) in different orders in the expansion in ^/l — z by a 
term (1 — a- 2)/-z, where the a^ ' are free parameters. Since the construction becomes singular for aj = —1, 
we exclude this value and use variations in the ranges a^ > and a>^' ' < —2. Note that for \a>^ ' '^j ^ oo the term 
{I — al ' ' z)/z becomes z-independent and the results for II^j^^ ' become independent of the a\ ' ' . We note that 
given the functions in Tab. E] it is straightforward to account for even higher terms in the expansion in the threshold 



limit for the construction of 11. 



(2,3) 
thr 



High-Energy Logarithms. The expansions of Il^^''^''^'{z) and G{z) in the high-energy limit \z\ — > oo read 



n(i'(z) = - 0.018998 log(-4z)- 0.075514 -0.056993 



ln(-4z) 



0.023628 



- 0.014248 



W(-4z) 



0.011874 



z 
ln(-4z) 



+ Oiz-'), 



n(2) (z) = (0.034829 - 0.0021109?!/) \n{-4:z) + (-0.050299 + 0.0029205n/) ln(-4z) + i7^^^ 

In^t—iz] \n(—4:z) h'-^^ 
0.18048 - 0.0063326nf) ^^ + (-0.59843 + 0.027441nf)— ^ + -^— 

z z z 

0.042745 - 0.0010554n/)^^-tM + (-0.10132 + 0.00580497i/)^^-tlfl 

\n(~4z) H^'^^ 
-0.48134 + 0.032065n/) ^ ' + -\- + 0{z-^) , 



n(3)(z) =(- 



0.063853 + 0.0077398n/ - 0.00023454 nj) ln''(-4z) 



0.21906 - 0.026441n/ + 0.0004867rii) In (-4z) 



-0.46209 + 0.10679n/ - 0.002183771^) ln(-4z) + W^ 

In^i—iz) 
-0.45120 + 0.035885n/ - 0.00070362n^) ^ - 



(3) 



(3.0848 - 0.26016n/ + 0.0045735n^) 



ln2(-4z) 



-6.6516 + 0.7823771/ - 0.014658771^)- 



•' z 

-0.10152 + 0.0060687n/ - 0.000087952n^) ^ 



(3) 



z z 

ln^(-4z) 



ln^(—4:z] In^f— 4zl 

0.57013 - 0.04485671/ + 0.0006743n^) — ^— — - + (0.17822 + 0.03852571/ - 0.0008839471^) — ^— — ^ 

Inf— 4z) H^^^ 
-8.8712 + 1.039371/ - 0.0260197i^)— -^-^ + ^- + 0{z-^) , 



G(z) 



log(-4z) 1 - log(-4z) 



2z 



4z2 



+ o(z-3). 



(17) 



The expansions for H^^^ and G are known from the exact expressions given in Eqs. (||) and (^, while the expansion 
for n^^) was taken from Ref. ||2^. Note that many orders in high-energy expansion are known for XI^^^ |g^, but we 
only consider in this work terms up to order 1/z^, since our analysis for II"' mainly serves as a testing ground for 
the application to II'^'. The expansion for If^^' was obtained in Rcfs. ||l^. At 0{a'^) the non-logarithmic coefficients 
Hq I are known analytically [g[ and read 



H, 



H 



(2) 


(2) 



0.73628 + 0.03764571/, 
0.30324-^0.02900271/. 



(18) 



At 0{a^) the non- logarithmic coefficients Hq ( have not been computed from Feynman diagrams in the Hterature 
before. As we show in Sees. VI] and VIII they can be determined from the reconstructed vacuum polarization function 

n(3). 



\n"{-4z) 


(1 - zr lG{z)r 


iln"(-4z),(n> 1) 


(i~zr-^[Giz)r 


iln(-4^) 


^G{z) 


^ln(-4z) 


'-^G{z) 


^ln^(-4^) 


^ lG{z)f 


jj\n\-4z) 


(i^ [G{z)]^ 


^ln^(-4^) 


{l-zf[G{z)f 



TABLE II: First column: Logarithmic terms that arise in the high-energy expansion of n"'"^^(z;) where 



oo. Second 



column: Corresponding functions used in the construction of IIj^^j ( 



(2.3)/ 



(2 3') 

To construct Hl^f {z) we have to find functions that can account for the different combinations of powers of ln(— 4^:) 
and of powers of 1/z that arise in the expansions of Eq. (jl^). A convenient choice is given in Tab. ||. Our guideline 
for including the factors of (1 — z)' is to ensure that the functions are constant or ~ y/1 — z in the threshold limit 
z ^ 1. This leads to 



Urj =B^'^'-±^il zfGizf + fsf^) + ^] i±^(f - z)Gizf 



t(2) 



-"0 



B. 



(2) , B.. 
30 ' 



^Vl-^)G(z) + i?f^i^^G(zf, 




+ B',t\l-z)'Gizf+\B,, 



(3) 



B 



(3)' 
41 



(f-z)G(z)^+ 5^0^ + ^^ + 



(3) 



B 



(3) 
51 



B, 



(3)~ 
52 



(f-z)G(z), 



(19) 



where the coefficients i?^ ' can be determined unambiguously from the conditions in Eqs. ( p^ ) and (^7|). In analogy 



t(2.3) 



^(2,3) 



to n^j^^ ' we have included modification factors (1 + 6^ ' "^ z)/z for the functions that are related to highest-power 



logarithmic terms at each order in the 1/z expansion. In n[^A^ we have a common modification factor for the 
functions related to the terms In (— 4z)/z and In (— 4z)/z^. For the parameters bl ' the choice 6^ = is excluded 
because in this case the construction becomes singular. For our analysis we adopt variations in the ranges \b[ 



(2,3) 



> 1. 



Using functions along the lines of Tab. g it is straightforward to account for even higher terms in the high-energy 



t(2,3) 



expansion in the for the construction of IIj^^j 

C0( 

the modification functions and then solves a set of linear equations 



Note that for the determination of the coefficient A, ' and B; ' one first fixes the constants a,- ' and 6, ' in 



Subtractions at q^ = 0. There are singularities ~ 1/z and ^ 1/z^ in II\y^j. (z) and IIj^^j (z) that arise in the limit 



z ^ 0. They are a consequence of the functions used to construct njj^j. (z) and H] 



(2,3), 



i,-jf (z). These singularities lead to 
unphysical behavior and need to be subtracted. For this task we define the function 



^^zero {^) — *-'( 



(2,3) 



s[ 



(2,3) 



S. 



(2,3) 



(20) 



After the coefficients Ai ' ' and B^ ' have been computed, the coefficients Sq { 2 arc determined such that 

<Y^(0) = 0. (21) 

(23).. . . . f2 3) 

Note that it is not mandatory to fix iSq in this way, and that our approach is independent of the choice for Sq . 
However, to satisfy Eq. (o) it is convenient for the purpose of presentation to impose the condition (pTj) and also the 
relation n^cg^^(O) = 0. 

V. DESIGNING H^og 

The terms Ilreg in Eq. ( |ll| ) have to account for the non-logarithmic conditions in the expansion at the threshold 
and at high energies, and for the coefficients that arise in the expansion around z = 0. We start by presenting the 
small-z expansion of II'^' and fl''^): 

n(2) = (0.719976 - 0.0296233n£)z + (0.698894 - 0.027 533 Ane)z^ + (0.637986 - 0.0240088n£)z3 

+ (0.584109 - 0.021162171^)2" + (0.539450 - 0.0189263^^)^^ + (0.502392 - 0.0171420n£)z'^ 
+ (0.471258 - 0.0156884nf)z^ + 0{z^) , 

n(3) = (10.6103 - 1.30278n£ + 0.0282783n^)z + (10.4187 - 1.12407nf + 0.022370671^)0^ + O^z^) . (22) 

The coefficients for 11'^^) were computed in Refs. R 0]. Recently the coefficients for 11'^^) have even been determined 
up to order z^'^ [pOl pl|. For II'^'^) the coefficient of order z was computed in Refs. [Ol O, and the coefhcient of 
order z^ was given in Ref. |16| . The coefhcients of order z" with ri > 3 have not yet been computed from Feynman 



diagrams. However, they can be determined from the reconstructed function H") as we show in Sees. VH and VHI. 



Designing Ilreg- We start exemplarily with the construction of Hrog. Close to threshold H^^^ exhibits the Coulomb 
singularity ~ l/Vl — z, see Eq. ([l5|). To avoid that the Fade approximant contains explicitly this singularity, we use 
two different methods: 

(i) We relate the Fade approximant P{lli) to /(z)Hrog, where f{z w 1) ~ vl--z. The coefficient of the Coulomb 
singularity is implemented through a condition on P(l). 

(ii) We use the relation 

and account for the Coulomb singularity by adding the function ^G{z) to Hj^j^. The Fade approximant P is 
not affected by the Coulomb singularity. 

The numerical differences that result from these two methods of implementing the Coulomb singularity constitute 
another tool for quantifying the uncertainties inherent to our approach. 

For method (i) the expression we use for the relation between the Fade approximant P{lj) and Hreg reads 

^(^) = 7TT^ l^rlli^^ - ^rlli-00)] , (24) 

where (^7^)2 ~ Vl ~ ^ for z -^■ 1. A similar relation was also used in Ref. B. Since the prefactor grows linearly with 
z, P(— 1) is a finite number. Some comments are in order concerning the term nrcg(— 00) that appears in Eq. (|24| ) and 
also in the analogous relations (|27| ) and ( p9[) that follow below. From the conditions H(0) = Hiog(O) = Hi-cg(O) = it 
is easy to see that 

P(0) = -H,.eg(-cx)). (25) 
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Thus in case that n].eg(— c)o) is known and taken as an input, Eq. ( |25| ) represents a condition that is imposed on the 
Pade approximant P. On the other hand, if nreg(— oo) is unknown or not taken as an input, it can be determined 



from Eq. (25) once the Pade approximant has been fixed from other conditions. We show in Sees. VII and VII] that 
this allows to determine the high energy constants H^ ' ' with small uncertainties. From Eqs. (p4h, (pq) and (ll^) the 
vacuum polarization function n"' is recovered from the relation 



n(2)(z) 



(1 + ^)^ 

l-UJ 



P{u) F(0) + ng(z) 



(26) 



From Eqs. (p3) it is now straightforward to determine the conditions on the Pade approximant P{u}) from the non- 
logarithmic constraints on n(z) in the threshold and the high-energy regions and from the coefficients in the expansion 
around z = 0. Additional constraints on P arise from the fact that in the limit |z| — > oo the first term on the RHS 
of Eq. ( P6[ ) can exhibit odd power terms ~ l/2;(2n+i)/2 -^^rj^j^ j^ _ 1^2,..., which do not exist in the high-energy 
expansion of the vacuum polarization function. It is reasonable to exclude such terms up to order l/z(2n-i-i)/2 ^j^g^ 
the information from the high-energy expansion up to order 1/z" is accounted for. For example, excluding terms 
~ 1/2'^/^ in Eq. ( P6|) leads to the constraint P(— 1) — 2P'(— 1) = 0, where P' refers to the derivative of P{u}) with 
respect to uo. Excluding also terms ~ 1/z^/^ leads to the condition 3P(— 1) — 9P"(— 1) -I- 2P"'(— 1) = 0. The various 
conditions on P lead to a complicated non-linear set of equations for the coefficients of the Pade approximant in 
Eqs. (|l2|), which we do not present explicitly here. These equations frequently have multiple solutions and are most 
conveniently tackled numerically. 

For method (ii), where the Coulomb singularity is treated in ^\' the relation between the Pade approximant P{uj) 
and Ivio],{z) reads 



P{lu) 



(1 



^?jAz) 



nL'l(-oo) 



(27) 



A similar relation was also used in Ref. g. Here, excluding terms of order 1/z'^/^ and 1/z^/^ for \z\ -^ 00 corresponds 
to the conditions P(-l) - -P'(-l) = and 6P(-1) - 6F"(-1) + P"'(-l) = 0, respectively. The vacuum polarization 
function is then recovered from the relation 



m^^(z) = {1+ujYP{lo) - P(0) + 



nS(^) 



(28) 



Designing Ilrog. 



The construction of Ilrog proceeds in a similar way. At 0{a'l) the vacuum polarization function has 



a Coulomb singularity ~ 1/(1 — z). For method (i) this singularity is incorporated in Ilrog and the relation between 
P{u}) and Ilrog reads 



P{lo) 



1 



n(3)(z)-n(3)(-oo) 



(29) 



where (1 — ujY 
conditions P'{— 
the relation 



(1 — z) for z ^ 1. Excluding terms of order 1/ z^/"^ and 1/z^/^ for \z\ -^ 00 corresponds to the 
= and 3P"(— 1) — P"'(— 1) = 0, respectively. The vacuum polarization function is recovered from 



n(3)(z) 



l-LO 



P{lu) - P(0) + 



n|S(.). 



(30) 



t(3) 



For method (n) we add the function 8C3[G(2)]V9 to U\'^^. Since 8C3[G'(z)]V9 -* 2 7r2C3/[9(l - z)] for z ^ 1 the 



t(3) 



(3) 



Coulomb singularity is therefore accounted for in ^\'- The relation between the Pade approximants and Ilreg has 



t(3) 



then the same form as Eq. (p4|) with Ilrcg replaced by Ilrcg. The relation for the n''^)(z) has the same form as Eq. (|2£ 



with n{^^ replaced by n[^^. The relations imposed on P to exclude terms of order 1/z'^/^ and Xj z^l'^ are then also 
P(-l) - 2P'(-1) = and 3P(-1) - 9P"(-1) -t- 2P"'(-1) = 0, respectively. 



11 

VI. DISCARDING UNPHYSICAL SOLUTIONS 

For the reconstructed functions II'^' and H^'^^ we have several types of variations that can be implemented into 
the construction and which we can use to quantify numericahy the uncertainty of the results. Apart from the two 
ways to account for the Coulomb singularity described as methods (i) and (ii) in the previous section, we have also 
implemented modification factors in Eqs. (Jl^) and ( |l9| ) that allow us to scan over a continuous set of functions within 
llj^ ' . Once the modification functions are fixed there are in general several possible choices one can use for the Fade 

approximants Pm.n with n + m being fixed by the number of conditions one imposes on Hreg ■ The resulting solutions 
for the Fade approximants can, however, have properties that lead to an unphysical and pathological behavior for II 
and R. Such solutions need to be discarded for a meaningful phenomcnological analysis 0] . 

An obvious restriction concerns solutions for Pm,nii^) that lead to poles in II'^''^-' in the complex w-plane inside 
the unit circle.^ These solutions are unacceptable and we discard them right away because such poles lead to an 
unphysical analytic structure. A more subtle situation arises for solutions with poles in n*-^'"^' in the upper complex 
w-half-plane that are outside the unit circle, but are either close to the unit circle or have a large residue. Although the 
analytic structure of such solutions is not a priori wrong, we still discard such solutions if they lead to an unphysical 
resonance-like structure in the cross section R. To have a quantitative criterion that can be implemented easily 
automatically we compute for every pole in Iliog in the upper complex w-half-plane the so called pole factor 

|Resn(c^poie)| ,^^. 

I^^polcl — 1 

where cjpoic is the location of the pole in the complex w-planc and Rcsn('j-'poic) the residue of II"''^^ at Wpoio. If l^poicl 
is close to unity or if the residue is large, the pole factor becomes big and a resonance-like structure can arise in 
R. We discard solutions when p > pQ. For our analysis we found that the choice po = 2.8 represents a reasonable 
restriction for 11^^^ while for n*^^^ we use po = 30. For the vacuum polarization function ll'^'^) a larger value for po is 
used since such poles arise predominantly in the threshold region close to a; = 1. Here Imp^'^^J is substantially larger 
than Im[n(^^] due to the bigger size of its Coulomb singularity, see Eq. (p^. Given the set of solutions for ll^^^ that 
pass the restrictions described above we can analyze how well these solutions reproduce other well-known properties 

ofn(2). 

VII. ANALYSIS FOR THE VACUUM POLARIZATION AT Oia^) 

The purpose of this section is two-fold. First we demonstrate the reliability of our approach for its application to 
n^'^^ by testing it with the rather well-known 0{a1) vacuum polarization function ll"^ and, second, we determine the 
previously unknown constant A' ^^^ that appears in the nonrelativistic expansion of 11^^^ close to the threshold, see 

Eq. ®- 

To demonstrate the reliability of our approach let us reconstruct 11^^^ using only information from the different 

expansions that is analogous to the available information in ll''^) . Thus we account for the expansions in the threshold 
region up to NNLO, in the high-energy region up to order 1/z^ and up to order z^ for the expansion around z = 0. 

(2) 

For the construction of Ilrog this entails that we account for the first two coefficients of the expansion around z = 0, 
the non-logarithmic term ex y/1 — z in the threshold region and the constraints that terms ~ 1/z^/^ and ~ Xj z^l"^ are 
absent for \z\ —^ oo. We do not implement the known constants Hq {, but we determine them from the reconstructed 
n"^ This amounts to 6 constraints on the Fade approximant for method (i), where the Coulomb singularity is 
accounted for in Hreg, and to 5 constraints on the Fade approximant for method (ii), where the Coulomb singularity 

(2) 

is accounted for in Xll'. Thus we have n + m ^ 5 for the Fade approximants Pm^n for method (i) and n + m = 4 for 
the Fade approximants for method (ii). 

Given the analytic form for the reconstructed 11'^' functions we can expand them in the threshold region, the 
high-energy limit and around q^ = 0. In Fig. |l| the results for the coefficients Cj, for k = 3,4,5,6,7, the high- 
energy constants Hq { and the threshold constant K^^^^ are displayed exemplarily for n/ = n^ + 1 = 4 relevant for 



For Fade approximants of the form Pk.o (Taylor-liko) such poles do not exist and none of the solutions is discarded. 
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■ 3.795 ± 0.079 



1: Results from the reconstructed II'^' function in approximation C for the coefficients Cg 4 5 g ^ 



that arise in the expansion 



r(2) 



FIG. 

around z — 0, for the coefficients Hq'^1 that occur in the non- logarithmic terms in the high-energy expansion \z\ -^ cxd and for 
the constant K^'^' that appears in the expansion at threshold around z = 1. The red solid lines represent the respective exact 
results known from computation of Feynman diagrams. The dashed blue lines represent the envelope of all results obtained 
from the reconstructed Tv ' functions. The individual error bars represent the range of vales obtained from the reconstructed 
n'^' functions using one particular Fade approximant Pm,n- The various types of Fade approximants that have been used are 
indicated in the upper left panel; the same order is used for all other panels. All results are for n/ — nt + 1 = 4: running flavors 
relevant for charm production. 



the production of charm quarks. The labels [m, n] which have been added to the upper left panel for C3 ' refer to 

(2) 
the Pade approximant used for the respective Ilreg functions, and their order is representative for all diagrams. The 

error bars represent the range of values covered by the variations of the modification factors as described in Sec. IV. 

The blue dashed lines indicate the range covered by all individual results and the red solid lines show the exact result 

obtained from Feynman diagrams. We see that for all cases the exact values are well within the range covered by the 

reconstructed 11^^^ functions. Particularly precise determination are obtained for Cg and the leading high-energy 

coefficient H^ ' . We also obtain a very precise determination of the the threshold constant /f^^' . 





approx. A 


approx. B 


approx. C 


approx. D 


approx. E 


exact 


^(20) 












2.49671 


C^^°' 












2.77702 


cf°) 


1.365 ±0.425 


1.609 ± 0.266 


1.611 ±0.048 






1.63882 


cf°) 


0.283 ± 0.799 


0.770 ± 0.441 


0.750 ±0.085 






0.79555 


^(20) 


-0.389 ±1.057 


0.271 ±0.521 


0.225 ±0.106 


0.278 ±0.001 




0.27814 


cr 


-0.744 ±1.213 


0.021 ±0.541 


-0.047 ±0.115 


0.007 ± 0.002 




0.0070080 


c|^°) 


-0.871 ± 1.296 


-0.054 ± 0.528 


-0.136 ±0.117 


-0.086 ± 0.003 


-0.08594 ± 0.00003 


-0.085963 


rr(2) 
-"0 


0.159 ±0.770 


-0.561 ±0.063 


-0.580 ±0.008 


-0.5854 ± 0.0004 


-0.5857 ±0.0001 


-0.58570 


J/f) 


0.007 ± 0.574 


0.338 ±0.871 


-0.132 ±0.071 


-0.180 ±0.008 


-0.185 ±0.004 


-0.18723 


kw 


-6.64 ±10.12 


3.933 ± 0.303 


3.795 ±0.079 


3.809 ± 0.032 


3.805 ± 0.020 





TABLE III: Results for the coefficients Cg 4 g g 7, Hq [ and K^^' from the reconstructed II''^' functions using various different 
types of approximations for II' '. Empty entries for coefBcients Cj. indicate that they are exact in that particular approxi- 
mation. The exact analytic form for K^'^' is unknown. All results are for rif — nt + 1 = A running flavors relevant for charm 
production. 



To demonstrate that our approach is systematic we need to show that the results become more accurate once more 
information is included for the reconstruction of H^^^. In Tab. Ill the results of Fig. |l| for C^ ' with k = 3,4, 5, 6, 7, 



Hq { and A'(^^ are displayed in the line labeled as "approximation C" . In comparison we also show the results when all 
the information ~ 1/z^ in the high-energy expansion is neglected for the reconstruction of H^^^ ("approximation B") 
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and when in addition to that also all NNLO threshold information is neglected ( "approximation A" ) . The results show 
that the properties of the vacuum polarization can be determined more accurately once more information is used for 
its reconstruction with our approach. Moreover, we find that the variations of the reconstructed vacuum polarization 
function due to the different choices for the modification factors and the Pade approximants represent a reliable tool 
to estimate the uncertainties. 



r>(20) _ f-<(20) 
0.05 



0.03 
0.01 
-0.01 
-0.03 
-0.05 
-0.07 



-0.09 




FIG. 2: Difference of the values for Cf, , determined from the reconstructed II"^ functions based on Taylor-like Pade 
approximants, and the known exact values Cj. ^,^acf ^^^ results are shown for approximation C (green triangles and light greed 
shaded area), D (red squared symbols and medium red shaded area) and E (blue diamonds and dark blue shaded area). The 
shaded areas represent the variation of the results due to the changes in the modification factors. 



^(20) 



It is an astounding and amusing fact our approach allows for a determination of the coefficients Cj, for large values 
of k with practically negligible uncertainties. This is demonstrated in Fig. y where the difference of the results we obtain 



for the coefficients (7^ and the exact values from Ref. pOl ^, Cj. 



Cl j,^^^j are shown up to fc = 30. The range 



of values between the green triangular-shaped symbols (green light shaded region) is obtained from approximation C 
using solutions with Taylor-like Pade approximants. For k > (10, 15, 20) the maximal relative discrepancy to the exact 
values is below (14, 3, 1)%. The discrepancies become even smaller when more of the coefficients in the expansion 
around z = of Eq. ( ^ ) are accounted for in the reconstruction of II'^^. Including the coefficients up to order z'^ 
(approximation D) we obtain the range of values between the red squared symbols (red shaded region). In this case 
we obtain for k > (10, 15, 20) a maximal relative discrepancy to the exact values of below (3, 0.8, 0.5)%. Including 
the coefficients up to order z^ (approximation E) we obtain the range of values between the blue diamond-shaped 
symbols (blue dark shaded region). Here, we obtain for fc > 10 a maximal relative discrepancy to the exact values of 

(2) n — 

below 0.2%. The results for the high-energy coefficients Hq { for approximations D and E are also shown in Tab. III. 



Again we find agreement with the exact results with decreasing uncertainties once more information is included for 
the reconstruction of II^^-' . Given the excellent quality of the results we consider our approach a reliable method to 
determine the previously unknown threshold constant K'''^\ As our final result for K'^'^^ wc adopt 



iC^^) ^ 3.81 ± 0.02. 



(32) 



To conclude this section let us analyze the 0{a1) corrections to the e+e cross section obtained from the recon- 



structed n(2). In Fig. I we have plotted l2-Kvlm[T\^'^\q^ 
pole mass scheme as a function of the quark velocity v 



+ zO)] for 71/ = 4 relevant for charm quark production in the 



yjl — \/z. We have included the factor of v to suppress the 
Coulomb 1/f -singularity that arises in the cross section for small values of v and to have a finite value in the limit 
iJ — > 0. In the left panel the result for approximation C is shown. The red band is the area covered by all solutions for 
n^^^ that pass the criteria discussed in Sec. VI and represents the uncertainty. The size of the uncertainty corresponds 



to the envelope of the individual error bars shown in Fig. ^ where approximation C has been used as well. For compar- 
ison we have also displayed the expansions in the threshold region for w — > (dotted lines) and in the high-energy limit 
for V ^> \ (dashed lines), where the short lines refer to leading order, the medium- length lines to next-to- leading order 
and the longest lines to next-to-next-to-leading order. The uncertainties are reduced substantially when additional 
coefficients for the expansion around z = are included for the reconstruction of n*-^' . This is demonstrated in the 
right panel, where the coefficients up to order z^ are included for the reconstruction of II"' . Again the red band 
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1271 V Im[n(2)] 

100 




127uvlm[n(2)] 

100 




FIG. 3: Results for 127r?;Iin[n"'((jr'^ + iO)] as a function of v for n/ — 4. The red bands represent the uncertainties. In the 
left panel the results are based on the reconstructed n'^' function incorporating the coefficients in the small-2 expansion up 
to order z (approximation C) and in the right panel the coefficients up to order z are accounted for. The dotted and dashed 
black lines show the expansions in the threshold and the high-energy region up to NNLO. See the text for more details. 



is the area covered by all solutions for IX^^^ that pass the criteria discussed in Sec. Vl. For method (i) to account 
the Coulomb singularity we found solutions based on the Fade approximants [9,0], [8,1], [7,2], [6,3], [5,4], [3,6], [1,8], 
and for method (ii) we found solutions based on the Fade approximants [8,0], [7,1], [6,2], [5,3], [4,4], [3,5], [1,7]. The 
width of the band is already smaller than the width of the solid lines used to draw the boundaries of the band. For 
V = (0.2, 0.4, 0.6, 0.8) the relative uncertainty is ±(0.09, 0.4, 2.0, 2.5)% and thus negligible for all conceivable practical 
applications. The approximation formulae for the 0{a'^Cp) and O^alCACp) contributions given in Ref. ||, g^l (using 
Eqs. (65) and (66) of Rcf. [D) together with the analytically known fermionic corrections agree within 1-2% with our 
result. We thus confirm the results for the cross section given in Refs. m, EG]. 



VIII. ANALYSIS FOR THE VACUUM POLARIZATION AT 0(af) 

For the reconstruction of H*^^^ we use all available information from the expansions in the threshold region, Eqs. (p^). 



the high-energy region, Eqs. (|7|) and around z = in Eqs. (22). For the construction of Ilrog we account for the 
first two coefficients in the expansion around z ~ 0, the non-logarithmic term oc \/l — z in the threshold limit and 
the two constraints from the absence of terms ~ l/z'^/^ and ~ 1/z^/^ for \z\ -^ oo. This amounts to 6 constraints 
on the Fade approximants for method (i), where the Coulomb singularity oc 1/(1 — z) is accounted for in Ilrcg; and 5 
constraints on the Fade approximants for method (ii), where this Coulomb singularity is accounted for in ^\'- Thus 
we have n + m = 5 for the Fade approximants Pm,n for method (i) and n + rn = 4 for the Fade approximants Pm.n 
for method (ii). 

In Fig. ^ the results for the coefficients Cf. for fc = 3, 4, 5, 6, 7, the high-energy constants H^ { and the threshold 



constant K^^^ are displayed for n/ 



nf 



1 = 4 relevant for charm quark production. The different labels [m, n] 



"(30") / . . (3) 

which have been added to the upper left panel for Cg refer to the Fade approximant used for the respective Ilrog 
function and their order is representative for all panels. The error bars represent the range of values covered by the 
variations of the modification factors as described in Sec. IV, and the blue dashed lines indicate the range covered 
by all individual results. We adopt this range as the uncertainty in our determination of these coefficients, and the 
results are summarized together with the corresponding results for n/ = n^ 4- 1 = 5 in Tab. IV. For the determination 



of the high-energy coefficients H^ and H{ we find uncertainties of about 1% and 10%, respectively. This compares 
well with the corresponding results for H^ ' and H{ ' we have obtained at 0{a^^) for approximation C, see Fig. |l| and 



Tab. |II|. For the coefficients Cj, with fc > 3 we find somewhat larger relative uncertainties than in for the C\. in 



approximation C. This is, however, not unexpected since the cancellations that arise when the pole mass results for 
these coefficients are transferred to the MS mass scheme are substantially larger at C'(aj). The result for K^'^^ has 
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FIG. 4: Results from the reconstructed H'"^^ frinction for the coefficients Cg 4 g g ^ that arise in the expansion around 2 = 0, for 

the coefficients Hq ( that occur in the non-logarithmic terms in the high-energy expansion \z\ — > cx) and for the constant K^^' 
that appears in the expansion at threshold around z — 1. The dashed blue lines represent the envelope of all results obtained 
from the reconstructed II'^- functions. The individual error bars represent the range of vales obtained from the reconstructed 
n'^' functions using one particular Fade approximant Pm.n- The various types of Fade approximants that have been used are 
indicated in the upper left panel; the same order is used for all other panels. All results are for n/ = n^ -|- 1 = 4 running flavors 
relevant for charm production. 



a particularly large error and can merely serve as a rough constraint on its true values. Concerning the precision in 
the determinations of K^^") , we believe that a substantial improvement can be achieved once the full set of NNNLO 
terms cx y/1 — z in the expansion for R at the threshold and the exact values for Cf. with A; > 3 become available. 





71/ = 4 


71/ = 5 


^(30) 


-5.6404 


-7.7624 


(^(30) 


-3.4937 


-2.6438 


^(30) 


-3.279 ± 0.573 


-1.457 ±0.579 


(^(30) 


-4.238 ±1.171 


-1.935 ±1.201 


(j(m 


-4.996 ± 1.666 


-2.507 ±1.732 


Mm 

<^6 


-5.280 ± 2.045 


-2.809 ±2.150 


^(30) 


-5.151 ±2.321 


-2.847 ± 2.467 


rr(3) 
^0 


-6.122 ±0.054 


-4.989 ± 0.053 


Hf> 


-3.885 ±0.417 


-3.180 ±0.405 


;^(3) 


-10.09 ±11.00 


-5.97 ±10.09 



TABLE IV: Summary of the results for the coefficients C3 4 5 g 7, Hq { and K''^^ obtained from the reconstructed II'^^ function 
for 71/ = Tif + 1 = 4 and ti/ = n^ + 1 = 5. The coefficients C} 2 ^^^ known exactly and shown for completeness. 



One of the most important applications of the coefficients Cn is the determination of the MS charm and bottom 
quark masses from moments M„ of the charm and bottom quark e^e~ cross section. For small values of n one way to 
compute the moments is using fixed-order perturbation theory as shown in Eq. dlQ). Using the results from Tab. IV 
we find for the fixed-order moments at 0{a'l) for charm quarks (ri/ = 4) 

Nh = (0.1348 ± 0.0044 ± 0.0005) x lO'^ , 
M4 = (0.153 ± 0.032 ± 0.002) x 10"^ , 
Ms = (0.199 ± 0.084 ± 0.008) x lO"'' 
Mq = (0.084 ± 0.144 ± 0.036) x 10"^ . 
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FIG. 5: Result for 127ri;IrQ[n'^'(g'^ + iO)] as a function of v for n/ = 4 using the currently available information for the 
reconstruction of n*-"^^ .The red band represent the uncertainty. The dotted and dashed black lines show the expansions in the 
threshold and the high-energy region up to next-to-next-to-leading order. See the text for details. 



Here we used 7nc{mc) ~ 1.27 GcV for the MS charm mass and aj^^ (1.27 GcV) = 0.387637 for the strong coupling as 
the input and four- loop rcnormalization group evolution. The first error arises from the variation of the rcnormalization 
scale between 1.27 and 3.81 GeV and the second error is due to the uncertainties in the 0{al) coefficients C^ shown 

5) with mbimb) = 4.17 GeV and ai"^=^^(4.17 GeV) = 0.224778 as the input 

Ms = (2.350 ± 0.017 ± 0.002) x 10"^ , 
Mi = (2.167 ± 0.045 ± 0.005) x 10"^ , 
M5 = (2.126 ± 0.091 ± 0.011) X 10"" , 



in Tab. IV. For bottom quarks (ny 
we find 



Me = (2.160 ± 0.148 ± 0.022) x 10 



-13 



(34) 



The first error arises from the variation of the renormalization scale between 2.085 and 8.34 GeV and the second error 
is due to uncertainties in the 0{al) coefficients Cj. shown in Tab. fv|. The results show that that the uncertainties 



in M3 4 5 caused by the errors in the coefficients Cg 45 we have obtained in this work are an order of magnitude smaller 
than the overall uncertainties of the moments at 0(0;^) due to variations of the renormalization scale. For physically 
relevant values of n they can be safely neglected. 

Finally, let us analyze the 0{al) corrections to the e+e" cross section obtained from the H^'^^ In Fig. we have 
plotted the function 127rulm[n'-^-' {q^ + iO)] for rif = 4 relevant for charm quark production in the pole mass scheme as 
a function of the quark velocity v = y/1 — 1/z. As for the analysis in Fig. H we have included the factor v to suppress 
the Goulomb singularity. The function still diverges logarithmically for v -^ because the 0{al) cross section has a 
singularity ~ ln(i')/w in the nonrelativistic limit. The red shaded band is the area covered by all solutions for II'^' 
that pass the criteria discussed in Sec. [V| and represents the uncertainty. The relative uncertainty is about 10% at 
V = 0.2 and 0.8 and should be acceptable for most applications where 0{a^) accuracy is required. For comparison 
we have also displayed the expansions in the threshold region for v ^ (dotted lines) at NLO (short line) and at 
NNLO (long line). Likewise the expansions in the high-energy limit for w — > 1 (dashed lines) are shown, where the 
short line refers to order l/z*^. the medium-length line to order 1/z and the longest lines to order 1/z^. We strongly 
emphasize the importance of incorporating the NNLO contributions in the expansion close to the threshold and the 
1/z^ terms at high energies for achieving our result. Once more information from the different kinematic regions 
becomes available, the uncertainties can be further reduced substantially. 
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IX. CONCLUSIONS 

111 this work we have determined the fuU mass and q^ dependence of the 0{a1) and 0{a'l) corrections to the heavy 
quark vacuum polarization function Ii{q^ , m^) and its contribution to the e+e^ total cross section. Our approach uses 
known results for the expansions of I{{q^,rin?) at high energies, in the threshold region and around g^ = 0, conformal 
mapping and the Fade approximation method. We have demonstrated for the vacuum polarization function at 0{a'j.) 
that the approach allows for reliable determinations of other properties of 11 with small uncertainties, and that the 
uncertainties of the results can be systematically reduced if more information from the three different kinematic 
regions is accounted for. Our results for the cross section at ©(a^) also confirm previous results by Chetyrkin, Klihn 
and Steinhauser from Refs. [|[ |2^. For the vacuum polarization function at 0{a1) we have determined the previously 
unknown non-logarithmic constant term that arises at NLO in the expansion close to the threshold. For the C(a^) 
corrections to the vacuum polarization function we determined the previously unknown coefficients in the expansion 
around q^ = Q beyond order q^ and the first two non-logarithmic coefficients in the high-energy expansion. The results 
for the coefficients in the expansion around q^ = allow for the determination of the moments M„ of the e+e~ cross 
section for n > 3 at C(aj). 
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Note added: After completion of this work K. Chetyrkin pointed out to us that analytic expressions for the constants 
H{ can be derived from results given in Ref. |g^. Evaluated numerically they give H{ ~ —4.33306 for nj = A and 



(3) ... . . h- tI 

HI = —3.53165 for Uj = 5 which is in agreement with the results we have presented in Tab. IV. 
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